Optical pumping in dense atomic media: Limitations due to reabsorption of 

spontaneously emitted photons 
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Resonant optical pumping in dense atomic media is discussed, where the absorption length is 
less than the smallest characteristic dimension of the sample. It is shown that reabsorption and 
multiple scattering of spontaneous photons (radiation trapping) can substantially slow down the 
rate of optical pumping. A very slow relaxation out of the target state of the pump process is 
then sufficient to make optical pumping impossible. As model systems an inhomogeneously and 
a radiatively broadened 3-level system resonantly driven with a strong broad-band pump field are 
considered. 



I. INTRODUCTION 



Optical pumping is an established technique in atomic and molecular physics to selectively populate or depopulate 
specific states or superpositions . It is based on the absorption of photons of a specific mode and subsequent 
spontaneous emission into many modes. The dissipative nature of the latter part makes it possible to transform mixed 
into pure atomic states. From this results the importance of optical pumping for state preparation in systems with a 
, thermal distribution of population and for laser cooling || . 

The maximum achievable rate of pumping is determined by the escape time of the emitted photons, which in 
optically thin media is given by the free-space radiative lifetime. When the medium becomes optically thick, however, 
i.e. when the absorption length becomes smaller than the smallest sample dimension, the escape time of photons can 
be substantially reduced. This phenomenon, known as radiation trapping [Q, is due to reabsorption and multiple 

■ scattering of spontaneously emitted photons and can drastically reduce the rate of optical pumping in dense media. 
These limitations could be of major importance in many different fields as for instance near-resonance linear and 
nonlinear optics in dense media |^,|6| or the realisation of Bose condensation by velocity selective coherent population 
trapping (VSCPT) §■ 

To describe the reabsorption and multiple scattering of photons we here utilize a recently developed approach to 
radiative interactions in dense atomic media Q . In this approach a nonlinear and nonlocal single-atom density matrix 
equation is derived which generalizes the linear theory of radiation trapping Q to the nonperturbative regime. As a 
model system a 3-level A configuration driven by a strong broad-band field is considered and the limits of (i) large 
inhomogeneous and (ii) purely radiative broadening are studied. 

Let us consider the A-type system shown in Fig. 1. A strong driving field with (complex) Rabi-frequency £l(t) 

■ couples the lower state |c) to the excited state \a), which spontaneously decays into \c) and \b). Since \b) is not 
coupled by the driving field, this results in optical pumping from \c) to \b). We also take into account a possible finite 
lifetime of the target state described by a population exchange between the lower states at rate 70. 




FIG. 1. Optical pumping in a A system. 



It was shown in M that the effect of the incoherent background radiation can be described by additional (nonlinear 
and nonlocal) pump and relaxation rates and level shifts in the single-atom density matrix equation. If we assume 
orthogonal dipole moments or sufficiently different frequencies of the two optical transitions, the level shifts are 
negligible. Also if the driving field is strong, the incoherent photons do not affect the pump transition a <-» c. Thus 
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we are left with a pump and decay rate T(t) on the a «-> b transition and the effective single- atom equations of motion 
read in a rotating frame: 

Paa = -(7 + 7' + r)/9 aa + Tp bb + i{fl*p ac ~ C.c), (1) 
Pec = l' Paa + loPbb - loPcc ~ «(^* Pac ~ C.c), (2) 
Pac = -(«A ac + r QC )/9 ac + ifl(p aa - Pec)- (3) 

A ac is the detuning of the drive field from resonance and r ac is the respective coherence decay rate. It should be 
noted, that T is a function of the density matrix elements of all other atoms, and hence the Eqs. are nonlinear 

and nonlocal. 

We are here interested in genuine optical pumping and therefore consider a broad-band pump || , i.e. f2(£) is assumed 
to have a vanishing mean value and Gaussian 5-like correlations (f2*(i)fi(t')) = R8(t—t'). Formally intergating Eq.(||), 
substituting the result back into Eqs.(|l]) and (||), and averaging over the Gaussian distribution of the pump field leads 
to the rate equations 

Paa = -(7 + 7' + F)p aa + T ' Pbb ~ R{Paa ~ Pec), (4) 
Pec = i Paa + loPbb - TOPcc + R(Paa ~ Pec)- (5) 



II. COLLECTIVE DECAY RATE 



We now have to determine the collective rate T. T is proportional to the spectrum of the incoherent field at the 
position r and the resonance frequency oj of the atom under consideration § 

r(u,t) = ^D(t f ,u ] t) = ^ rdr((E-(fo,t)E + (fo,t + r)))e^. (6) 

tl It J — oo 

Here E ± are the positive and negative frequency parts of the field operators, p is the dipole matrix element of the 
atomic transition, and ((AB)) = (AB) — (A)(B). D(oj) can be obtained by summing the spontaneous emission 
contributions of all atoms propagated through the medium S 

D(l,l) = JJ d3d4D rct (l,3) (D rct (l,4))*n s (3,4). (7) 

Here D ret (l,2) is the retarded propagator of the electric field inside the medium, which obeys a Dyson-equation in 
self-consistent Hartree approximation: 

D™ t {l,2) = Dl ct {\,2)- JJ rf3d4L>f ) ct (l,3)n rot (3,4)D rct (4,2). (8) 

In Eqs.(Q) and (||) the numbers 1,2... stand for {fx, t\}, {ra, ia} . . ., and the intergrations extend over time from 
— oo to +oo and over the whole sample volume. Dg ct is the free-space retarded propagator of the electric field. For 
simplicity we here have disregarded polarisation. We also have introduced the atomic source correlation 

n s (fx,t i; f 2 ,t 2 ) = ^^((J j (t 1 )a j {t a )))6(r 1 -f j )5(ft-? j ) (9) 
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and the atomic response function 
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n ret (ri,*i;r 2 , ta) = ^ &(ti - t a ) ^( [a}(h), aj(h)] ) 6(n - r-) 6(r 2 - f 3 ), (10) 

3 

where dj — \b)jj{a\ is the spin-flip operator of the jth atom and O is the Heaviside step function. In terms of 
the rr's the dipole operator of the jth atom reads dj — p((Jj + crt). The names reflect the physical meaning of the 
quantities (|^, |l0| ) . The Fourier-transform of n s is proportional to the spontaneous emission spectrum of the atoms 
and that of n r5F gives the susceptibility of the medium. Eqs.(^|) and (0) represent a nonperturbative summation of the 
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spontaneous radiation contributions of all atoms propagated through the medium. It assumes a Gaussian statistics, 
which is however a good approximation for the background radiation. 

The Dyson-equation (§h was solved in with some approximations in a macroscopic (continuum) limit where 
II(fi, fi; r*2, £2) = / d 3 rP(r, t±, £2) 8(r\ — r) 5{f t 2 — r)- This yielded for the collective decay rate 

where r = l^o — r| is the distance bewteeen the source and the probe atom. The probability that a photon reaches 
the probe atom is determined by the absorption coefficient 

4{r, w,t) = ^- Re [P ret (f, u; tj\ . (12) 
One can easily calculate the atomic source and response functions for the A-system of Fig. 1. 



pLajt) ~ Pj b (t) 

h 2 r ab + i(w - U} J ab ) 



P rct (r 3 ,u,t) = J ,aaw / Ffcb 7 x , (13) 



p S(f t) = V^ pL(*)r ot , 

J & a (r Qb ) 2 + (c-^' 

where N is the density of atoms, tu ab is the resonance frequency of the jth atom, T ab the coherence decay rate of the 
corresponding transition and the overbar denotes averaging over a possible inhomogeneous distribution of frequencies. 

At this points we shall distinguish two limiting cases. We first consider the limit of large Doppler-broadening and 
secondly the case of purely radiative broadening. 



III. INHOMOGENEOUSLY BROADENED SYSTEM 



The approach of Q is based on the Markov approximation of a spectrally broad incoherent radiation. This ap- 
proximation is justified for example in an inhomogeneously broadened system. We therefore discuss first the case of 
large Doppler-broadening. If we are interested in the population dynamics on a time scale slow compared to velocity 

changing collissions, we may set plAt) = pjj,^(t) = P^^{rj,t) and thus have the same population dynamics in all 
velocity classes. Since Y depends on the populations of all atoms, Eqs.([|) and (|j) are nonlocal. In the case of a 
constant density of atoms and a homogeneous pump field, Y and hence all density matrix elements will be approx- 
imately homogeneous. We therefore make a simplifying approximation and disregard the space dependence. The 
volume integral is then carried out by placing the probe atom in the center of the sample. This yields for a Gaussian 
Doppler-distribution of width 3> 7 
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Paa(t) 



Pbb(t) - Paa(t) 



l-exp(-ff(t) e - A2 / 2A -) 



(15) 



where A = to — uj^ b is the detuning from the atomic resonance at rest, and H (t) = K [pbb(t) ~ p aa (t)]- K — g NX 2 d c g 
with g = j /\^2n Ad characterizes the number of atoms within one relevant velocity class in a volume given by the 
wavelength squared and the effective escape distance d e g. In deriving (|l5|) we have used the relation between the 
free-space radiative decay rate 7 and the dipole moment p: p 2 — Sirheoc 3 "/ /uj 3 ||l0| . d c g corresponds for a long 
cylindrical slab to the cylinder radius; for a thin disk to its thickness and for a sphere to its radius. 
Averaging over the inhomogeneous velocity distribution of the atoms eventually yields 

30 -1 
dLu^L e- A2 / 2A -lV,t) 

Paa(t) 1 



r(t)= Y(u,t) 



7 



Pbb{t) - Paa(t) y/lT 



dye 



-v 



1 — exp 



(-H(t)e 



(16) 



In Fig. 2a we have shown the population in the target state \b) as function of time starting from equal populations 
of levels |c) and |6) at t = 0. We here have assumed that the target state is stable, i.e. 70 = 0. One recognizes that 
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optical pumping is considerably slowed down already for values of K on the order of 10, which usually corresponds 
to much less than one atom per A 3 . The slow-down of pumping is further illustrated in Fig. 2b, where the effective 
pump rate defined as 



dt 



In [p aa + pc 



(17) 



is plotted normalized to the value in an optically thin medium (r° = 7/2). One can see that the optical pump rate 
approaches a constant asymptotic value, which for K 3> 1 and large pump rates R is given by 



2K{n\nK) 



7 

1/2^2 



(18) 



p bb (t) 



1 






0.9 










1 i 




O.S 








0.7 




'i 






i 

/ 


1 






/ 


5 


0.6 




/ 


10 






/ 


50 


0.5 









los r p -1 



-1.5 
-2 
-2.5 







K= 1 


\ 

V 

\ 


5 


\ 


10 


50 



100 200 300 400 500 100 200 300 400 500 

ft yt 

FIG. 2. (a) Time evolution of population of level \b) for 70 = 0, R/'y = 10 and 7Y7 = 1 for different density parameters 
K — g N\ 2 d c s , g = 7/v / 27rZ\ c . (b) Corresponding effective rate of optical pumping 



Since we have assumed in the plots of Fig. 2 an infinitely long-lived target state (70 = 0), all populations eventually 
ends up in \b). However if 70 is nonzero and in particular if it becomes comparable to the asymptotic rate r^f, the 
steady-state populations of all states equalize. In this case optical pumping is less and less efficient and becomes 
eventually impossible. This is illustrated in Fig. 3, where the stationary population in state |6) is shown as a function 
of the density parameter K for different values of 70. 
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FIG. 3. Stationary population in level |6) for R/'y = 10, V/7 = 1 and different values of 70 as function of density parameter 
K. 



IV. RADIATIVELY BROADENED SYSTEM 



We now discuss the case of a radiatively broadened system. In analogy to the case of inhomogeneous broadening, 
we find for the spectral distribution 
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Paa{t) 



Pbb{t) ~ Paa{t) 



1 - exp -H{t) 



lab^- ab 



(19) 
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where A = u -co ab , T ab = 7 ab + r, and -f ab = (7 + 7' + R + 7o)/2, and = K [p bb (t) - p aa (t)]- Here K = gN\ 2 d eS 
with g = r )/2n^f ab . As opposed to the corresponding relation in the inhomogeneous case, Eq.([l9|) determines the 
collective decay rate only implicitly, and T needs to be calculated self-consistently. For small atomic densities or 
Paa ~ Pbb the exponential function in Eq.(|l9|) can be expanded into a power series. The first nonvanishing term 
found from this has the same spectral shape than the single-atom response function. In such a case the Markov 
approximation used in Q is no longer valid and the approach is quantitatively incorrect. We shall nevertheless use it 
and discuss the range of validity afterwards. 

We find that in the case of radiative broadening the rate of optical pumping decreases exponentially with the density 
parameter as opposed to [NX 2 d c s]~ 1 in the inhomogeneous case. For sufficiently large pump rates R and stable target 
state (70 = 0) the asymptotic rate of optical pumping is here 

T™ = lexp{-K}. (20) 

Physically this is due to the fact that here the incoherent photons are in resonance with all atoms, which drastically 
increases the scattering probability. As a consequence much smaller decay rates 70 out of the target state are sufficient 
to make optical pumping impossible. This is illustrated in Fig. 4, where we have plotted the stationary population in 
state \b) as function of the density parameter K = N\ 2 d c ff /2ir for different values of 70. 
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FIG. 4. Same as Fig. 3 for radiatively broadened system; Ko = NX 2 d c s /2w, R/"f = 10, 7V7 = 1 

In order to check the validity of the Markov approximation, we have shown in Fig. 6 the stationary normalized 
spectral distribution T(lu)/T for K = 1, 10 and 100 and 70/7 = 10~ 4 . Also plotted is the atomic absorption spectrum 
for Kq — 1 (solid line). 
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FIG. 5. Spectral distribution of incoherent background radiation for R/j = 10, 7' = 7, 70/7 = 10 -4 , and Ko = 1 (dotted), 
Ko — 10 (dashed) and Ko ~ 100 (dashed-dotted). Also shown is the normalized absorption spectrum for Ko = 1. 

One recognizes that spectrum of the background radiation has only a slightly larger width than the atomic response 
for Kq = 1. In this case the Markov approximation is not valid. The situation however improves when the density is 
increased. Thus Fig. 4 has only qualitative character for lower densities. 





5 



V. SUMMARY 



We have shown that resonant optical pumping in a dense atomic medium is substantially different from optical 
puming in dilute systems. When the absorption length of spontaneously emitted photons process becomes less than 
the minimum escape distance, these photons are trapped inside the medium and cause repumping of population. This 
leads to a considerable slow-down of the transfer rate and can make optical pumping impossible if the target state 
of the pump process has a finite lifetime. The effect is much less pronounced in inhomogeneously broadened systems 
due to the reduction of the spectral density of background photons. 

These results may have some important consequences. It is practically impossile to use resonant optical pumping 
in media with NX 3 ~ 1 . This sets strong limits to the possibility to prepare pure states or coherent superpositions in 
systems with initial thermal occupation of states, such as Hyperfine ground levels of alkali at room temperature. Even 
though the above analysis did not take into account quantum properties of the atoms and considers only resonant 
pumping, the results indicate, that it may be very difficult to achieve Bose Condensation via VSCPT in optical lattices 
jllj . Also the present results show that electromagnetically induced transparency (EIT) |Q in dense media cannot be 
understood as the result of optical puming into a dark state. Essential for EIT in dense media is an entirely coherent 
evolution jl3| via stimulated adiabatic Raman passage |I3]. Some of these aspects will be discussed in more detail 
elsewhere. 
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